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Abstract. We consider front solutions of the Swift-Hohenberg equation d^u = —{\ + d^)u + eu — u. These 
are traveling waves which leave in their wake a periodic pattern in the laboratory frame. Using renormalization 
techniques and a decomposition into Bloch waves, we show the non-linear stability of these solutions. It turns 
out that this problem is closely related to the question of stability of the trivial solution for the model problem 
d^u{x, t) = d^u(x, t) + {I + tanh(a; — ct))u(x, t) + u{x, t)^ withp > 3. In particular, we show that the instability 
of the perturbation ahead of the front is entirely compensated by a diffusive stabilization which sets in once the 
perturbation has hit the bulk behind the front. 
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1. Statement of the problem 

We consider the Swift-Hohenberg equation 

d^u = -{I + dl)\ + e\ - u\ (1.1) 

with u{x, t) G R, .X G R, t > and < e ^ 1 a small bifurcation parameter. It has been shown 
some time ago that a 2-parameter family of (small) spatially periodic solutions exists which are 
independent of t. These solutions correspond to a periodic pattern which exists in the laboratory 
frame. These solutions are of the form 

Ug,aix) = cos((l + eq)x + a) + 0{e^) , 

which bifurcate from the solution u = 0. Here, 

It is furthermore well-known and proved in [CE90a] that these solutions are marginally stable 
for 4|g| < ^, the so-called Eckhaus stability range ([Eck65]), and that the spectrum of the 
linearization about these solutions is all of R~. Finally, after a long time it was shown in 
[Schn96] that these solutions are also non-linearly stable, and this proof was the presented in a 
slightly different form in [EWW97]. 

In another direction, in earlier work of [CE86] and [EW91] traveling wave solutions of 
a special kind leaving a fixed pattern in the laboratory space were shown to exist, and their 
linear stability was studied in [CE87]. Our present paper is concerned with a first proof of the 
non-linear stability of these traveling solutions. 

We first describe the traveling solutions. One way to view them is to write 

ne2Z+l 

where A^^ ~ A^ as defined above and A^ _^ = A^ ^, with x the complex conjugate of x. 

Here, the A^ ^ are in fact 0(5'"'), and furthermore ^ extends to an analytic function. The 
modulated front solutions are then of the form 

U{x,t) = F^ q^^{x - Ct,x) , 

with 

Fc,ai^x)^^, 5^ W^^^^^{C)e^-(('+^'^^-+^K (1.2) 

n£2Z+l 

Note that these are not classical traveling waves of the form u{x — ct), and note furthermore 
that i^c,q,o periodic in its second argument (with period 27r/(l -|- eq)). The modulated front 
solutions satisfy [CE86, EW91], when c> 0: 

lim W.^^^^ii) = A^^^ , lim = • 
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These modulated front solutions are constructed with the help of a center manifold reduction, 
where all W^c,g,n determined by the central modes VFc,q,±i- reduced four-dimensional 

system for q ±i = q ±i(0 there is a heteroclinic connection lying in the intersection of 
a four-dimensional stable manifold of the origin and a two-dimensional unstable manifold of an 
equilibrium corresponding to ^. Since this is a very robust situation these solutions can be 
constructed by some perturbation analysis from the ones for q = 0. For small e and q = the 
solution Q J of the amplitude equation on the center manifold is close to the real- valued front 
solution i(0 — ^B{6^) = eB{Q of the equation 

4dlB + c^d^B + B- 3B\B\'^ = , 

connecting o,i = at C = +oo with o,i — C = —oo. The constant is given 

by Cb = s~^c = 0(1). Our paper deals with the question: Under which conditions does the 
solution of (1.1) with initial data ^^(a;, x) + v{x) converge to F^^^^^{x — ct, x) as t oo? 

We will show our results for the case q = and a = only, to keep the notation on 
a reasonable level. The extension to arbitrary a is trivial by translating the origin, while the 
extension to arbitrary q satisfying 4\q\ < ^ necessitates some notational work and leads to 
bounds which depend on q. Thus, we will write the periodic solution as 

U^{x) = Acosx + O{e^) , (1.3) 

with A = 2e, and the modulated front (moving with speed c = 0(e)) as 

nG2Z+l 

We describe next the nature of the stability problem. Consider an initial condition Uq{x) = 
F^{x, x) + Vq{x), and let u{x, t) denote the solution of (1.1) with that initial condition. Since 
solves (1.1), we find for the evolution of v{x, t) = u{x, t) — F^{x — ct, x): 

d^v{x,t) = {Lv){x,t) - 3F^{x - ct,xfv{x,t) - 3F^{x - ct,x)v{x,tf - v{x,t)^ . (1.4) 

Here, L = —(1 -|- d^)^ + e^. We define the translation operator r^^ by {T^^f){x) = f{x — ct, x), 
so that (1.4) can be written as 

d,v = Lv- 3{T,,Ffv - 3{t,,FJv' - . (1.5) 

Introduce now K^^ (the difference between the modulated front and the periodic solution) by 

Kjx) = (t,,F,) (x) - (x) = F^ix - ct, x) - t4 (x) . (1.6) 

Note that K^^{x) vanishes as x — > — oo, and approaches U^{x) as a; — > oo. With these notations 
we can rewrite (1.5) as 

d,v ^ Lv- 3U^v - 6t/, K^,v - 3Kl,v - 3U, - - 3Ky 
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where 

Mv = Lv- 3U^v , 

The variables with index i vanish with some exponential rate for fixed x G R in the laboratory 
frame. They will be seen to be exponentially "irrelevant" in terms of a renormalization group 
analysis. In order to explain this renormalization problem, we will study, in the next section the 
model problem 

dfu{x, t) = d^u{x, t) + a{x — ct)u{x, t) + u{x, t)^ , 

with a(^) = ^(1 + tanh^), and p > 3. This problem is nice in its own right. The similitude 
will come from the correspondence of M with d^, and of M-^v with the term a{x — ct)u{x, t). 
Indeed: 

• the first term will be seen to be diffusive in the laboratory frame, 

• the second term will be seen to be irrelevant in the laboratory frame, but the first together 
with the second term will be exponentially damping in a suitable space of exponentially 
decaying functions in a frame moving with a speed close to c. 

As in previous work [Sa77, BK94, Ga94, EW94] our analysis will be based on an interplay of 
estimates obtained in these two topologies. 

Our main results are stated in Theorem 4. 1 for the simplified problem and in Theorem 7. 1 
for the Swift-Hohenberg problem. We not only show convergence to the front, but give also 
precise first order estimates in both cases. As far as possible, the treatment of the two problems 
is done in analogous fashion, so that the reader who has followed the proof of the simplified 
problem should have no difficulty in reading the proof for the full, more complicated, problem. 

Remark. An ideal treatment of this problem would necessitate a norm in a frame moving with 
the same speed as the front. Such a space is needed to study the stability of so-called critical 
fronts (moving at the minimal possible speed where they are linearly stable). Achieving this aim 
seems to be a necessary step in solving the long-standing problem of "front selection" [DL83], 
in a case where the maximum principle [AW78] is not available. 

Remark. The method also applies to more complicated systems, like hydrodynamic stability 
problems. A typical example are the fronts connecting the Taylor vortices with the Couette flow 
in the Taylor-Couette problem. These fronts have been constructed in [HS99]. The stability of 
the spatially periodic Taylor vortices has been shown in [Schn98]. 

Notation. Throughout this paper many different constants are denoted with the same symbol C. 

Acknowledgements. Guido Schneider would like to thank for the kind hospitality at the Physics 
Department of the University of Geneva. This work is partially supported by the Fonds National 
Suisse. The work of Guido Schneider is partially supported by the Deutsche Forschungsge- 
meinschaft DFG under the grant Mi45 9/2-3. 
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Part I. A simplified problem 

2. The model equation 

Let 

a(0 = i(l+tanhO. (2.1) 

We want to study the equation 

d^u{x,t) = dl.u{x,t) -\- a{x — ct)u{x,t) -\- u{x,tY , (2-2) 

with c > and p>3. For notational simpUcity we assume p e N. 

To understand the dynamics of (2.2) it might be useful to consider the following simplified 
problem 

dfv{x,t) = dl.v{x,t) + 'd{x — ct)v{x,t) , (2.3) 

where ^{z) — 1 when ^ > and d{z) = when ^ < 0. If we go to the moving frame — x — ct 
and let w(^, t) = v{x — ct, t), then the equation for w becomes 

dMC, t) = dlw{^, t) + cd^wiC, t) + mnC, t) ■ (2.4) 

For a; > 0, we have '&{x) = I and hence the corresponding characteristic polynomial for (2.4) 
(in momentum space) is 

-k^ + ick + l , 

while for a; < 0, we have '&{x) — with its corresponding polynomial 

—k^ + ick . 

Thus, we expect the solution to be exponentially unstable ahead of the front, i.e., for a; > 0, 
and diffusively stable behind the front. If we consider an initial condition f()(0 localized 
near ^ = > 0, and of amplitude A, then we expect the amplitude to grow like e^A until 
t — — ^q/c, when this perturbation "hits" the back of the front (in the moving frame), or, in 
other words, when the back of the front hits the perturbation (in the laboratory frame). Thus, 
the perturbation does not grow larger than Ae^''^^. We use this in the following way. Assume 
that the amplitude at ^ > is bounded by Ae~^^. Then, ignoring diffusion, we find that the 
contribution to the amplitude at the origin at time t = ^q/c is bounded by 

f^dC^e^^^-^^^/^ 

Clearly, if /?c > 1, the initial perturbations are sufficiently small for the total effect at the origin 
(in the moving frame) to be small. 

Once this has happened, a second epoch starts where the perturbation is behind the front. 
Then, due to the diffusive behavior, the amplitude will go down as 

C 

{t-K + iy^ ■ 

These considerations will be used in the choice of topology below. 
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2.1. Function spaces and Fourier transform 

We start the precise analysis and will work in Fourier space and revert to the a; -variables only at 
the end of the discussion. We define the Fourier transform by 



i^f){k) = ^ Jdxf{x)e- 



■ikx 



Notation. If / denotes a function, then / is defined by / = JF/, and if A is an operator, then 
A is defined by ^ = !FAT~^. We also use the notation f * g for the convolution product 
f9={f* 9){k) — J^^f{k — £)g{£)- Finally, denotes the conjugate of translation: 

(fjm = e-<^~m , (2.5) 
so that the Fourier transform of (T^f) (x) = f{x — Q is 

The relation ([Ta97]) 

k-d^,{j-f){k)^{-i)fnd>^f)ik) 

motivates the introduction of the following norms: We fix a small 5 > and define 

/ 2 \l/2 

Jdk\di{k'm)n . (2.6) 



The dual norm to this is 



/ 2 \ 

= [ ^52(^+.) Jdx\dif{x)\'x'^\ . (2.7) 



Parseval's inequality immediately leads to: 

UWhi, = II-^/IIh^.^, 

and, for some constant C independent of 1 > 5 > 0, 

WfgWnl, < C\\fh2 \\g^ 

I,o I,d Z,d 

Finally, we shall also need the inequality 

II/*^IIh^.^ < II/IIc^II^IIh^.^, (2-9) 



(2.8) 
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where 

2 

\\f\\ci, = J2^'^^p\dim\. (2.10) 

This follows from 

II/*^IIh^.^ = II/-^IIh2 < Il/llc2 II^IIh^ = Il/llc2 II^IIh^.^ , 

2 2,o D,d 2,0 D,o *^^2 

where the inequality above is a direct consequence of the definition of H2' . 

Notation. In the sequel, we will always write || • || instead of || • ||g2,5. Thus this is our default 
norm. ^ 
We define the map by 

(VV/3,ct/)(0 = f{C + &t)e^^, (2.11) 

where /? e (0, P^) and c e (0, c) will be fixed later. The Fourier conjugate of this operator then 
satisfies 

0^0, J) (k) ^ {:FW^,,t^-'f) (k) = e'^'+^^^^'fik + i/3) , (2.12) 
as one sees from the following equalities: 

M^^pjm = /dee-^"^(W^,,,/)(0 

= y"dCe-*('^+*^)«-'^*)/(0 
= 27re*('=+*^)^7(^ + iP) ■ 

This calculation also shows that if f{$.)e^*^ e Hl gfor f e C^g, then Wp^^f extends to an 
analytic function m {0 > ImA; > — and {y^p^^tf){' ~ ^/^) ^ for all (3 e [0, (3^). 

Remark. Since the norms for different 5 are equivalent, all theorems throughout this paper can 
also be formulated in a version with 5=1. 
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3. The linear simplified problem 

In this section we study the Unearization of equation (2.2): 

dtU{x,t) ^ dlU{x,t) + a{x-ct)U{x,t) . (3.1) 
The function a is given as 

a(e) = i(l+tanhe), (3.2) 

but our methods will work for many other functions. The crucial property we need is the 
existence of a > such that a(^)e~^^ satisfies 

U^a{Oe-^%. < C, (3.3) 

1,0 

for all (3 E (0, (3^). For the case of (3.2) we can take (3^ = 2. The Fourier transform a of a 
is therefore a tempered distribution which is the boundary value of a function (again called a) 
which is analytic in the strip {z \ > Imz > P^}. Furthermore, there is a K such that, for all 

S e (0, 1], 

||a||c2 < 1 + K5, (3.4) 

b,6 

since 

sup \a{x)\ < 1 . (3.5) 

The bound (3.5) will be tacitly used later. 

The next proposition describes how solutions of (3. 1) tend to as t — > oo. We write (x) 
for U{x, t) and use similar notation for other functions of space and time. 

Proposition 3.1. Assume that there are a P and a c G (0, c) such that jS^ — jdc + \ = —2'y < 0. 
Then there exists a S G (0, 1] such that the following holds. Assume that G H2 ^ and that 

Wq{C) = (>V^,o^o)(6 = ^o(?)e^^ e nj^s- (These conditions are independent of 5 > 0.) 
Then the solution U^{x) = U{x,t) of (3.1) with initial data Uq exists for allt>0 and with 

'^{k) = the rescaled solution V{k, t) = U{kt~^^^, t) satisfies 

\\V, - UomWnls < ^^^^WUohr ■ (3.6) 

The function — V^j^^ct^t satisfies 

\\W,\\^2,s < Ce-'^'/'\\Wj^2.s . (3.7) 
The constant C does not depend on Uq. 



Remark. Note that it is optimal to choose c arbitrarily close to c. 
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Proof. First of all, we rewrite the equation (3.1) for in terms of and Wf The equation for 

= >V^,,,C/, is 

d,Wi^, t) = ajWi^, t) + ic- 2/3)d^Wi^, t) 

+ a(e - (c - c)t)Wi^, t) + {(3^ - f3c)W{^, t) . 

Taking Fourier transforms, we then find, omitting the argument k and using the notation of (2.5): 

d,U, = -k% + (7;,a) * [/, , ^ ^ ^ (3.9) 

d,W, = {P^ -Pc-k^ + ik{c - 2(3))W, + (T(,_e),a) * W, . (3.10) 

It is at this point that the simultaneous choice of two representations for the solution and their 
associated topologies is crucial. 

We first show that converges to 0, i.e., we show (3.7). We find from (2.9): 

\\{ffa) * /II < ||a(- - Ollc,.^ • 11/11 = llallc,^^ • 11/11 • (3.11) 
Therefore, (3.4) implies 

||(T(,_g),a)*;^^,|| < {l + K6)\\W,\\ , 
and we get from (3.10) the bound 

for a constant independent of 5 E (0, 1]. The term K^S comes from the derivatives in the 
norm || • |U2,5. We choose 5 > so small that 

(3^ ^ I3c+1 + {K + K,)5 < -37/2 . 
Integrating over t we get from the choice of /?, 5, and c: 

\\W^\\ < e-^^*/^\\WQ\\ . (3.12) 

Thus, we have shown Eq.(3.7). 

Next, we study U. From (2.12) and deforming the contour of integration, we get 

((T^a) * f){k) = J d£e-^^(^-^)a(A; - £)/{£) 

= J dle-'^^^-^^a{k - (){yVp,J){£ - iP)e-'^^' 

= J d£e-'^^^-^-'^^~a{k - e - ip){W^^^J){e)e-'^^+ 

= e-^^^-^*) j dt e-'^^^-^^a{k - i - i(3) [Wp^gJ) {tfe'' 



I3)ct 



(3.13) 
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Let h{k) = e-'''"'a{k - 1(3) and g{k) = e-^'='^*(>V^^g,C/,) (A;) = e-'''^'Wt{k). Then (3.13) 
implies _ _ 

From this we conclude that 

\\{f,ia)*U,\\ =e-^^'^-'^'\\h*~g\\ 

< e-^^'^-^^i/ill ll^ll (3.14) 

< C{l + tc)\l + tcfe-l^^''-^^'\\Wt\\ . 

On the other hand, from (3.7) we know that \\Wf\\ stays bounded (it actually decays exponen- 
tially), and thus the evolution equation for is of the form 

dMk) = -k%{k) + h{k,t){l + tcf{l + tcfe-f^^^-^^\ 

with t)|| uniformly bounded in t. Since, by construction, c < c, we conclude that (3.6) 
holds, using well-known arguments which will be made explicit in the proof of Theorem 4.1. 
The proof of Proposition 3. 1 is complete. □ 



4. The renormalization approach for the simplified problem 

We consider now the non-linear problem (2.2) and its related version for = Wp g^Uf 
TyV^ ct^t Fourier space. It takes the form 



Qt^t = -k {T^th) * fi, + -u*P 



, n (4.1) 

d,w, = {p^ -pc-k^ + ik{c - 2/3))*, + (T(,_g),a) * *, + fi^^ 'Uw^. 

Let Mp be the operator of multiplication: {M^f){x) — e^^f{x). Choose the constants c, 
and P such that they satisfy as before 

> -27 = /^^ - /3c + 1 , 

and fix them henceforth. Our main result for the simplified problem is: 

Theorem 4.1. There are positive constants R, C and 5 E (0, 1] such that the following holds: 
Assume IIwqIIh^ + II-^a^oIIh^ — -R- Then the solution of (2.2) with initial condition 

2,6 ^ 2,6 

Uq converges to a Gaussian in the sense that there is a constant — A^{uq) such that with 

~ 2 

'il^{k) — e~'^ therescaled solution v{k, t) — u{kt~^^^, t) satisfies 

CR 

\\vt-Aj^s < (4-2) 
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Furthermore, 

\\wt\\^2,s = \\y'W^,,M^2,s < CRe--^' . 

We shall use the renormalization technique of [BK92] to show that and behave (as 
t — > oo) essentially in the same way as their linear counterparts C/^ and from the previous 
section. This technique consists, see [CEE92], in pushing forward the solution for some time 
and then rescaling it. This process makes the effective non-linearity smaller at each step, so that 
in the end the convergence properties of the linearized problem are obtained. 

We fix < cT < 1 and introduce: 

(£/)(x) = />x). (4.3) 
This is again a linear change of coordinates in function space. Note that 

= y"dx7>x-x')^(x') 
= a{{l~f)*{l~g)){H). 

Furthermore, 

(£(f^a))(x) = e<'^-~a{a>c) = (T;^(£a)) (x) , 

and therefore we have 

£((f^a) * /) = a{f^^L~a) * (£/) . (4.5) 

We next define 

w^^^{k) = e-'^""'""(£"w)(x,cT-2V) = e-^'""'"^w((j"x, (j-^V) , 
so that this corresponds to an additional rescaling of the time axis. Note that 

^A"") = e-^'^""^'^(a"x,a-^V^) = <_i,i(^-^x) . 

We also let = t'^a. From (4.4), (4.5), and = a~^'^df we find easily that (4.1) transforms 
to the system (omitting the argument x): 

= ((/32 - /3c + 7)^-^- + zx(c - 2/3)a--)<,^ (4.7) 
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We see that under these rescalings the coefficients of the non-linear terms go to as n — > oo. 
We will now put this observation into more mathematical form. 

The equation (4.1) is of the form dfX^ = L(^X^^ + 3\r(X^), where L contains the linear 
parts with the exception of those depending on and 'N denotes the other terms. We can write 
the solution as 



Jtn 



Going to the rescaled variables ^, and taking = cr ^^'^ ^^andt = a ^"r, we can express 
this (for the u) as follows. The equation (4.6) leads to 



(4.8) 



Similarly, we rewrite (4.7) as 
where ^ is defined, cf. (4.7), by 

{Gn,rf)H = ((/3'-/3c+7)<T-^--x2+ix(c-2/?)<7-'^)/(x)+<7-'^((f(,_,)^_„,aJ*/)(x) . 

The solution of the linear evolution equation d^f^ ^ = ^f^ ^ is nothing but (3.10) in a new 

coordinate system. We write the solution as ^ = ^ ^,f^ Then, in analogy to (4.8) we 
get 

<,rM = iSn,r,a^<,a^) (><) + ^'^^""'^ jjr (^S^^^^^, «^7^) * " (4-9) 



Remark. The proof of Theorem 4.1 is divided into several steps: In Lemma 4.2 below, we 

improve first the inequalities for the exponentially damped part in scaled variables. Then in 
Lemma 4.4 a priori estimates for the solutions of (4.8) and (4.9) are established. With these 
a priori bounds we show Proposition 4.5. From these results. Theorem 4.1 will follow rather 
simply by a contraction argument. 
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4.1. The scaled linear problem 

Here, we derive the essential bounds on the influence of the term a{x — ct)u{x, t) in the equation 
for under the scaUngs introduced above. Note first that, from definition (2.6) and (4.3), we 
have 



From this we conclude immediately that for < u < 1 : 



ll^/l < a-'/^\\f\\ and \\t-'f\\ < a-'/'\\f\\ . (4.10) 
We next bound S^^^,. Recall that we are assuming — /3c + 1 = —27 < 0. 
Lemma 4.2. For all e' e (0, 1) there exists aC^, > such that fori > t > t' >0 one has 

\\Sn,r,r'f\\ < C,,^-^ V^-""^---')/^/!! , (4.11) 

for all n e N. 

Proof. We consider the equation d^f^ = G^.^f^, whose solution is = S'^^^,/^,: 

dj, = A J, + (7-"(f(,_g),_„,aJ * , (4.12) 
where is the operator of multiplication by 

A^(x) = -Pc + 7)(7-2^ - + zx(c - 2/3)(7-" . 
The variation of constant formula yields 



We now introduce the norm 



and its dual 



We use 



= ^5^^ Ux\x^f{x)f . 
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and 



2 , N -2n 

a T 



for every x > 0, where the g are constants independent of a depending only on x and 5. 
They have the property that Uni^_^Q C^g = I for fixed x- We choose x = 7/4 and find 



-2n, /, 



since 



II (T^^a J = a-"||aJ--C)^-VllHO,, 

< ^-1|an(--C)|U|l-^-V~|lH0 = ^-"l|anlk«ll/llH^.- 



Using II a |Uo = 1 and applying Gronwall's inequality to e ^c+57/4)ct "^y f || -2,5 we get 

b "0 

or equivalently, 

ll/rllH^'^ < ^,/4,Jl/r'llH^.^e(^'-^^+5^/^+^V4,.)<^-^"(--') . (4.13) 
' 

We choose 5 e (0, 1] so small that C'-y/4,5 < 7/4 + 1. This proves the assertion of Lemma 4.2 
for the Hq' norm. 

We next use the regularizing character of — to prove the bound in 62''^. Let q'(x) = x. 

Then 

IIq/JIh^.. < Ce(^'-^'^+^^/^)-""(^--')||g7.'llH^.^ 


+ C rdse(^'-^^+3'''/2V"'"(^-^) sup je-'^'^^-^^xj a'^" ||a„||co||/JU2,4 . 
Jt' >ceR " ° 

Using the estimate (4.13) for ||/s|U2,6 we get 

||g/,||g2.. < C||(l + |Q|)/,,||^2,.e(^'-^^+i+37/2)a--(.-/) j^^^^i^ _ ^/)i/2^ _ 
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2 / X , ' 



To bound the second power of q, choose e' e (0, 1). Then 

Jt' xER 
""b 

where ||5'|||^o,e' = sup^^g \g{x)\ + sup^^R \{J^~^{mg){x)\ with rh{k) = |1 + k^]'^ Clearly, 
||o„|Lo,£' is finite and using the estimate (4.13) to bound ||/s|U2,5, we get 

\\fUh2.s < C||(l + g')/,,|U2,.e(^'-^'^+i+3^/')"""("-^'V-^''^max(l, (r - r')^'""'^/') • 

"o 

Combining these estimates completes the proof of Lemma 4.2. □ 

Remark. It is easy to see that additionally the following holds: For all e', a E (0, 1) there exists 
a C^, ^ > such that for 1 > r > r' > one has 



-2n , /, 



ll'^n,.,.'/llH^.^ < C^',a^-' "e-^<^" )/2(r - rO"||(l + | • lY'-^'flk^.s , 

' ' 

for all n e Z. 

4.2. An a priori bound on the non-linear problem 

We now state and prove a priori bounds on the solution of (4.8) and (4.9). Finally these solutions 
will be controlled by proving inequalities for the elements of the following sequences. 

Definition 4.3. For all n, we define 

Pn = \\<,l\\ = IK ill . 

Moreover, we define 

K = sup ||w„^^|| and R"^ = sup |K^^|| . (4.15) 

re[a2,l] Te[o-2,l] 



Lemma 4.4. For alln eN there is a constant ry^ > such that the following holds: If p^_i, 
Pn-i, and a > are smaller than rj^, the solutions of (4.8) and (4.9) exist for all r e [a^, 1]. 
Moreover, we have the estimates 

K < Ca-"^pl_, + Ce-^"""< + Ca''^P-^\Rlf , (4.16) 
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and 

K < C(T-^/^-"'"p;^_i + Ca''^p-^-''\R:;^f-'R^ , (4.17) 

with a constant C independent of a and n. 

Remark. There is no need for a detailed expression for rj — rj^ since the existence of the 
solutions is guaranteed if we can show < oo and i?^ < oo. With (4.16) and (4.17) we have 
detailed control of these quantities in terms of the norm of the initial conditions and a. 

Proof. We start with (4.9). We bound the first term of (4.9) by using a variant of (4.11): First 
note that ('§n,^,^2<,^2)(x) = (^(4-i,Ta-M<-i,i)) (^)- Therefore, 

\\^i^n-l,ra-\l^n-l,l)\\ = ^"^^^ll-^n-l.ra-^l^-l,!) II 

c/T ' -2"/ 2, . l4.ioj 

< Ca-^/V-^ "e-^" ^/'|K_i,i|| . 

Therefore, we get for the first term in (4.9) a bound 

C(j-^/V-^'"p;^_i . (4.19) 
For the second term in (4.9), we get a bound 

/ -i\ / / — 2ti ^ / 2-. , 

C^n(p-3) / ^^f^-en^-^a -a ) /2 ^^uy-l 

< C(7"^p-3-^')(j^"(i?;;)p-ii?;^ < Ca''^p-^-''\Rif-^R';; . 

We next consider (4.8). The first term is bounded by 

||x^e-'(---')<_i,i(ax)|| 

< ||x^e-'(---')||.. ||x^^,_i,i(ax)|| (4.20) 

using (4.10). Using (3.13) and (4.5), the second term can be rewritten as 



= a-'^ rdr'e-'(---')((f,,_.„,,a)*C^,_.„,,)(^"^) 
= cT"^" / dr'e-^'" T (c-c)g->. (r-T ) 



—2n I —2n I 
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Using this identity, we get from the techniques leading to (3.14): 



2n I —In I 

n,T' 



(4.21) 



_2(n— 1) — 7T, 

< (7(j-6ng-(/3(c-c)+7)a ^u; <^ (^g-(/3(c-c)+7)<T ^ 

For the last term in (4.8) we get a bound 

C^-(p-3) [\r\RlY < Ca''^P-^\Riy . (4.22) 

The proof of Lemma 4.4 now follows by applying the contraction mapping principle to (4.8) 
and (4.9). For p^_i, Pn-i u > sufficiently small the Lipschitz constant on the right hand 

side of (4.8) and (4.9) in C([(T^, 1], 62'^) is smaller than 1. An application of a classical fixed 
point argument completes the proof of Lemma 4.4. □ 



4.3. The iteration process 

We next decompose the solution ^ for r = 1 into a Gaussian part and a remainder. Let 

2 

V'(x) = e~'^ and write 

where r„(0) = 0, and the amplitude is in R. We also define n : ^ R by 

n/ = /L=o- (4-23) 
Then (4.8) can be decomposed accordingly and takes the form 

fjx) = e-'(i--')f„_i(ax) (4.24) 
+ £dr' e-^(^--') (a-"(f,^_„^,a J * + a^^^-'^ u^^^^^j (x) 

+ e-'(i-^')A„_i^(c7x) - Aj{x) . (4.25) 
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Then we define p!^ = \\r^\\ and so < C{\A^\+ plJ. Our main estimate is now 

Proposition 4.5. There is a constant C > such that for a > sufficiently small the solution 
u of (2.2) satisfies for all n e N; 

l^n - ^n-il < Ce-^'-'^K + Ca-^'-'HKf , (4.26) 

Pn < p;_l/2+ce-^'^""< + cc7-(^-3)(i^-)^ 

< Ce-^-""p-_, + Ca^^^-'-^'\Kr-'K ■ (4.27) 



Proof. We begin by bounding the difference — A^_^ using (4.24). Observe that since we 
work in H2' , we have 

|ff/| < Cll/ll , (4.28) 

with C independent of 5. Thus, it suffices to bound the norm of the integral in (4.24). The first 
term in (4.24) is the one containing the translated term a„ and was already bounded in (4.21) 
while the second was bounded in (4.22). Combining these bounds with (4.28), we find (4.26). 

We next bound f„ in terms of f^_i, using (4.25). The first term is the one where the 
projection is crucial: For cr > sufficiently small, / e ^2'^ with /(O) = one has 



X I— > e 



/>x)|| < II/II/2. 



(4.29) 



~ 2 (5 

Indeed, writing out the definition (2.6) of H2' , one gets for the term with j 



y"dxe-2"'(i-'"'^|/((7x)|2 = (7-iy"d((7x)e-2"'(i-'^')(c7x)2 



/(CTX) - /(O) 



(7X 



Clearly, a bound of the type of (4.29) follows for this term by the assumptions on /. The 
derivatives are handled similarly, except that there is no need to divide and multiply by powers 

of (Tx since each derivative produces a factor cr. 

We now bound the other terms in (4.25). The first term is bounded using (4.29) and yields 
a bound (in ^i^'^) of 



Pn-l/2 • 

The second and third terms have been bounded in (4.21) and (4.22): 

Finally, the last term in (4.25) can be written as 

^ A,_i(e-'^i-')e-'-' - e-) + {A^_, - AJe 



(4.30) 



(4.31) 
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~ 2 (5 

The first expression vanishes and we get a bound (in H2' ): 

ll^nll < Ce-^"""i?;;; + Ca''^p-^\RiY . (4.32) 

Collecting the bounds (4.30)-(4.32), the assertion (4.27) for follows. Finally, the bounds on 
follow as those in Lemma 4.4. The proof of Proposition 4.5 is complete. □ 



Proof of Theorem 4.1. The proof is an induction argument, using repeatedly the above estimates. 
Again we write C for (positive) constants which can be chosen independent of a and n. Assume 
that R — sup^^jvj < oo exists. From Lemma 4.4 we observe for u > sufficiently small 



— 5/2— ne w 



" " 1 - C(7^(P-l-=')i?P 



K < 



Rl 



(4.33) 



1 -Ca'^iP-^^RP-^ 



Pn-l ' 



with a constant C which can be chosen independent of R. Using Proposition 4.5 we find 

Pn < c(Ki+p;), 

Pn < Ce p„_i + Ca ^P >a I p^_^ . 

Therefore, we can choose u > so small that for n > 3: (recall p > 3 and p e N) 

\K-K-x\ < Pn-i/io + a--'(|A_j + p;_j, 

Pn < 3p;_i/4 + p^i/10 + a--3K_i|, 

Pn < P^-l/lO. 

Thus, the sequence of converges geometrically to a finite limit A^. Furthermore, we find 
that lim^^^ pj^ = 0, and lim^^^ = 0. Since the quantities \AJ, p'^, p^ increase only 
for at most three steps the term CR^~^ in (4.33) stays less than 1/2 if we choose p[, 
p^ = 0(0""^), for an m > sufficiently large. We then deduce from (4.33) the existence of a 
finite constant R — sup^^p^ R^. Finally, the scaling of ^ implies the exponential decay of 
w^. The proof of Theorem 4. 1 is complete. □ 
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Part II. The Swift-Hohenberg equation 



5. Bloch waves 

Since the problem we consider takes place in a setting with a periodic background provided by 
the stationary solution of the Swift-Hohenberg, it is natural to work with the Bloch representation 

of the functions. For additional informations see [RS72]. 

The starting point of Bloch wave analysis in case of a 27r-periodic underlying pattern is 
the following relation 



u{x) 



J-i/2 ~Ti J-m 



(5.1) 



where we define 



{Tu){e,x) = u{e,x) = ^e^"^w(n + £) . (5.2) 

nGZ 



The operator T will play a role analogous to that played by the Fourier transform for the 
simplified problem of Part I. We will use analogous notation: 

Notation. If / denotes a function, then / is defined by / = Tf, and if ^ is an operator, then A 

is defined by ^ = T AT~\ 
Note that 



[ dx\u{x)\'^ = 2tt de dx\u{e,x)\^ . (5.3) 

Jr J- 1/2 JO 

This is easily seen from Parseval's identity: 



( dx\u{x)\^ = In [ dA;|«(A;)|^ 

R Jr 

= 271 V / de\u{n + e)\' 



= 2n 


ri/2 

/ de 

1-1/2 


E 

neZ 


= In 


n/2 

/ de 


r2TT 

i 




n/2 





The sum and the integral can be interchanged in (5.1) due to Fubini's theorem when u is in the 
Schwartz space S. 
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We shall use frequently the following fundamental properties (which follow at once from 
(5.2)): 

u{e,x) = u{e,x + 2TT) , (5.4) 

u{£, x) — u{—i, x) for real-valued u . 

Multiplication in position space corresponds to a modified convolution operation for the Bloch- 
functions: 

n/2 

{irv){e,x) = / de'u{e-e',x)v{e',x) = {u®v){e,x) . 

.1-1/2 

This follows from (5.4) and the identities: 



{irv){l,x) = l'dku{e + m-k)v{k)e''^'' 
r\/2 

= de' u{e + m - e' - n)v{e' + n)e' 

■^"1/2 m,n€Z 

Recalling the norm 

1/2 



{m—n)x^inx 



nl, - [Y^'^"^^'^ [dx\d-f{x)\'xA 



we now introduce 



V,tlo -^1/2 Jo 

We get from Parseval's equality 

for some C independent of 5 e (0, 1). Similarly, in analogy to (2.8), we also have 

||'in;|L2,5 = ||'u®{i|U2,5 < C||'u|L2,« 11-011-2,5, (5.5) 

2 2 2 2 

\\(uv){- -ip,-)\\^2^,s = ||(it ®i))(- - •)||g2,5 < C||it||g2,5 ||i)(- - •) 11^2,5 . (5.6) 
Finally, suppose / is a function in ^ (see (2.10) for the definition): Then, 

11/^11 = \\f®v\\^2,S < C\\f\y \\V\\^2,S , (5.7) 

2 2 2 

\\{f®v)i--t(3,-)\\^l,s < C||/||c2^ ||5(.-z/3,-)||hM. (5.8) 

Thus, apart from notational differences, we can work in the Bloch spaces with much the same 
bounds as in the spaces used for the model problem of the previous sections. 



1/2 



|2 
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6. The linearized problem 

We discuss here again tlie behavior of the linearized problem as in Section 3, but now for the 
Swift-Hohenberg equation. The discussion will again be split in an aspect behind the front and 
one ahead of the front. In Section 3, the behavior of the problem in the bulk behind the traveling 
front was diffusive by construction, and the only difficulty was to understand the role of the 
decay of a to (as e"''''^') as a; — > — oo. For the problem of the Swift-Hohenberg equation, 
the situation is similar, leading again to diffusive behavior. However, this observation is not 
obvious. Therefore, the first problem consists in showing the diffusive behavior. In order to 
obtain optimal results for the analysis ahead of the front, i.e., for the variable in the weighted 
representation, we use our approximate knowledge of the shape of the front. 

6.1. The unweighted representation 

In analogy with the simplified example, the linearized problem would be now 

^^v = Mv + MiV , (6.1) 

where M and Mj have been defined in Eqs.(1.7) and (1.8). By the analysis for the model 
problem we expect that the term M{u will be irrelevant for the dynamics in the bulk with some 
exponential rate. Therefore, it will be considered in the sequel together with the non-linear 
terms. As a consequence, the linear equation dominating the behavior behind the front is given 
by 

d^v = Mv . (6.2) 

We recall those features of the proof of diffusive stability of [Schn96, Schn98] which are relevant 
to the study of (6.2). 

In order to do this, we need to localize the spectrum of M. Since this is well-documented, 
we just summarize the results. As the linearized problem has periodic coefficients, the operator 
M = TMl ^ equals a direct integral J®di M^, where each acts on the subspace with fixed 

^ 2 (5 

quasi-momentum i in H2' . The eigenfunctions of are given by Bloch waves of the form 
e^^^Wi ^ with 27r -periodic Wi ^. The index n G N counts various eigenvalues for fixed £. For 
each £ e R (or rather in the Brillouin zone [— ^ , ^]) they are solutions of the eigenvalue equation 

{Mgw^){x) = -{1 + {ie + d^fywi{x) + s^wg{x) -3U^{x)w^{x) = f^^w^ix) . 
The spectrum takes the familiar form of a curve /ij {£) with an expansion 

l^,{£) = -c/ + 0(£3), 

and Cj > and the remainder of the spectrum negative and bounded away from 0. The 

eigenfunction associated with /ij(0) is d^U^ (x), reflecting the translation invariance of the 
original problem (1.1). There is an £q > such that for fixed £ e (— ^0) eigenfunction 
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(Pi{x) = I {x) of the main branch /ij (£) is well defined (and a continuation of d^U^ {x)) as £ 
is varied away from 0. Corresponding to this we define the central projections -Pc(^) by 

where (-, ■) is the scalar product in L^([0, Itt]) and (^^ the associated eigenfunction of the 
adjoint problem. We will need a smooth version of the projection in H2' . We fix once and for 
all a non-negative smooth cutoff function x with support in [— £q/2, £q/2] which equals 1 on 
[—£0/4, ^0/4]- Then we define the operators E^. and by: 

It will be useful to define auxiliary "mode filters" E^ and E^ by 

El{i) = x(^/2)Pe(^) , = l(^) - X(2^)n(^) • 

These definitions are made in such a way that 

which will be used to replace the (missing) projection property of E^. and E^. 

We next extend the definitions (4.3) of Section 4 to the Bloch spaces. To avoid cumbersome 
notation, we shall use mostly the same symbols as in that section. Thus, with cr < 1 as before, 
we let now 

(£'u)(x, x) = u{ax^x) . 

Note that here, and elsewhere, the scaling does not act on the x variable, only on the quasi- 
momentum x. The novelty of renormalization in Bloch space here is that since the integration 
region over the i variable is finite it will change with the scaling. Therefore, we introduce (for 
fixed 5 > 0), 

^a,p = I ll^ll^^,, < 00 '} (6-3) 



where 



\u 



, ,A /•1/(2t) r27r ^ ' , , 

E / dx6'^-+^^\dp:^u{£,x)\\l+fr 

n,n-0-^-V(2a) Jo 



For technical reasons we introduced a weight in the Bloch variable i. We will always write /C^ 
instead of /C^ ^. Note that T, as defined in (5.2) is an isomorphism between the space ^ and 
the space /Cj by (5.3) and the definition (6.3). 

Consider again the eigenfunctions (fi^ix). The function 

v,{e,x) = e^^^'^'ip,{x), 
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solves the equation 

Because of the nature of the spectrum ji^ {£), this solution satisfies 

Using this observation and the fact that the E'^-part is exponentially damped, the result will be 
Proposition 6.1. The solution of the problem (6.2) with initial data Vq satisfies: 

\\{£,x)^V,{£t-'/\x)-e-'^^''p^{0)V,{0,x)\\^^^^ < ^||i>ollH^.^ , (6-4) 
for a constant C > and allt > 1. Moreover, there is a constant 7_ > such that 

\\i£,x) ^ (4 V,)i£t-'/\x)\\,c,,^^ < Ce-^-'\\Vo^s , (6.5) 

for all t> 1. 

6.2. The weighted representation 

The weighted representation will be obtained by translating the effect of the transformation 
£i defined in (2.1 1) to the language of the Bloch waves. In accordance with our notational 
conventions, we set ^ 

and we get now, in analogy to (2.12), 

(W^,ct/)(A^) = e''^'+'f'^*f{£ + ip,x + ct). 
The equation (6.1), expressed in terms of g^v, then takes the form 

dt (>V^,ct^) = Kct (>V^,ct^) + ^i,^,ct K,ct^) , (6.6) 

with 

{Mf,^J){£,x) = {L,pf){£,x)-2>U^{x)f{£,x) + c{i{£ + if3) + d,)f{£,x), 

(^i,/3,ct/)(^,^) = -6^. {x){k^^®j){£,x) - 3ik^, ®k^, ®mx) . 

Some explanations are in order: L^^ is the operator — ( 1 + {d^ + ii — + e'^. The functions 
are just multiplications in the Bloch representation because they are periodic. More precisely. 
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one has {£, x) = C/^ {x)S{i) in the sense of distributions. The functions K^^ are derived from 
K^^ of Eq.(l .6) and are seen to be given by 

K,,{e,x) = (rKj(£,x) = e-'''''F^{e,x-ct,x)-u,ix)Sie), 

where the Bloch transform is taken in the first (non-periodic) variable of F^. 

In order to obtain optimal results for the analysis ahead of the front, i.e., for the variable 
in the weighted representation, we the recall some facts from the construction [CE86, EW91] 

of the fronts. 

For small £ > the bifurcating solutions u of the Swift-Hohenberg equation can be 
approximated by 

'4'{x,t,e) = sA{ex, eH)e^'' + C.C. , 
up to an error 0{e^), where A satisfies the Ginzburg-Landau equation 

dj^A = 4dxA + A- 3A\A\'^ , 

with X e R, T > and A{X, T) E C. See [CE90b, vH91, KSM92, Schn94]. This equation 
possesses a real-valued front Af{X, T) = B{X — c^T), where B{^) satisfies the ordinary 
differential equation 

AB" + CbB' + B- 3B\B\^ = . 

For |c^| > 4 the real-valued fronts of this equation are monotonic. These fronts and the trivial 
solution A = can be stabilized by introducing a weight e^^ * satisfying the stability condition 

^a(cb, /34 ) = 4/?j - + 1< , 

see [BK92]. 

Remark. Since B{^) converges at a faster rate to 1 /V3 for ^ — > — oo than to for ^ — > oo there 
will be no additional restriction such as (3.3) on . 

Remark. Our result will be optimal in the sense that each modulated front which corresponds 
to a front of the associated amplitude equation satisfying ^»^(c^,/34) < is stable. The 
connection between the quantities of the Ginzburg-Landau equation and the associated Swift- 
Hohenberg equation is as follows. We have c = ec^ + 0{s^), and P = eP^ + 0{s^). 

In order to prove this remark we write the modulated front F^ as defined in (1.2) as a sum 
of the Ginzburg-Landau part and a remainder 

F^(e,x) = 2eB{eOcos{x)+e^F^{^,x), 

where Fj. satisfies 

sup||F,(- + y,-)||c2 < C, 

for a constant C independent of £ G (0, 1) and 6 e (0, 1). Then we consider (6.6) which we 
write without decomposition as 

d^W = {L^f,W) - 3{7^) ® (r^) ®W + ciiie + iP) + djW . (6.7) 
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In order to control these solutions we use that the linearized system (6.6) evolves in such a 
way that during times of order 0(1/ e^) it can be approximated by the associated linearized 
Ginzburg-Landau equation 

d^A = 4{dj, - fA + csid^ -p^)A + A- B\2A + A) . (6.8) 

Theorem 6.2. For all Cq > 0, and > there exist positive constants Eq, C^, C2, and Tq such 

thatforalle e (0,£q] the following is true: For all initial conditions Wq with \\Wq\\^2,6 < CqE 

2 

there are a solution of (6.7) and a solution A^ of (6.8) with ||^q|U2,5 < Cj such that the 
function A^ approximates in the sense that 

\\W,- eT{A^,,_^^{x)e^^ + c.c.)\\^2,s < C^e^ 

for all t G {tq/s^, (tq + rj/e^]. Here T again denotes the map of Eq.(5.2) from a function f 
of X to its Bloch representation /(£, x). 

Proof. The proof of this is very similar to the case of the (non-linear) Swift-Hohenberg equation 
which was discussed in the literature [CE90b, vH91, KSM92, Schn94]. Our (linear) problem is 
in fact easier and the proof is left to the reader. □ 

For the system (6.8) we have the estimate [BK92] 

with lim^^o g^ic^, f3^,S) = g^ic^, f3^ ). The deviation of ^?^(cg, (3^ , 6) from ^^(0^, ) 
comes again from the derivatives of B and from the polynomial weight in the norm g. As a 
consequence of this estimate and of Theorem 6.2 we conclude that 

\\W^\\^2,5 < Ce^('^'^'^'^)(*-*')||W^^,||g2,. , (6.9) 

for a constant C and a coefficient g = g{c, e, 5). We can (and will) choose this constant g in 
such a way that (for £ — > 0): 

g{c,P,e,6) = e\g^{c^,PA,5) + o{l)). (6.10) 

We define g{c, P, e) = lim^^Q g{c, P, e, 5). 

Remark. The choice of a sufficiently small S > and £ > will allow us to prove the 
stability of all fronts which are predicted to be stable by the associated amplitude equation since 

lim^^^s)^Qe~^g{c,P,e,S) = ^^(cs,/?^)- 
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In the following we consider a modulated front with velocity c and a given (sufficiently 
small) bifurcation parameter £ > for which there are a /3 and a c G (0, c) which satisfy: 

^(c,/?,£) = -27<0. (6.11) 



Proposition 6.3. Suppose that the above stability condition (6.11) is satisfied. Then there is a 
5 e (0, 1] such that: There isaC < oo for which the functions = g^V^^ obey the bounds 

\\Wt\U2,6 < Ce-3^(*-")/2||t?JU2,5 . (6.12) 
"2 2 

As in the previous sections this result will have to be improved for the non-linear problem. 
Therefore, we skip at this point the proof, and will only deal with the improved version later. 

Thus, the linear problems (6.2) and (6.6) are the analogs of (3.9) and (3.10) and can be 
studied pretty much as in the case of the simplified problem, yielding inequalities similar to 
(3.6) and (3.7). 



7. The renormalization process for the full problem 

We assume throughout this section that the stability condition (6. 1 1) is satisfied. We prove here 
our main 

Theorem 7.1. There are a 5 > and positive constants R and C such that the following holds: 
Assume \\vq\\^2 +\\ MqVq \\^2 < R and denote by the solution of (1 .4) with initial condition 

2,5 2,5 

Vq. Let ^(£) = exp(— Cj£^). There is a constant = A^{vq) such that the rescaled solution 

vl{£, x) = Vf{it~^^^, x) satisfies 

CR 



Furthermore, 



- = l|W^,cAk.,^ < CRe--^'. (7.2) 



Remarks. 

• The inequality (7.1) really says that the difference 

v,{£t-'/\x)-A,e-'^^''d,U,{x) 

is small, where is the periodic solution (see Eq.(1.3)) of the Swift-Hohenberg equation. 
Expressed in the laboratory frame, this means that an initial perturbation Vq{x) will go to 
like 

/~~7r~ — 
Vt{x) ^ A^{vQ)J—exp{^)d^U^{x), 
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when t ^ oo, uniformly for x G R. See [Schn96]. In particular, this means that near the 
extrema of the convergence is faster than 0{t~^^^) since at those points d^U^ vanishes. 

• The inequality (7.2) gives some more precise bound on the growth of a perturbation ahead of 
the front, because it says that this perturbation decays exponentially in the weighted norm. 
More explicitly, we have at least a bound 

\vt{x + ct)\ < Ce^^-^'* , 

with 7' slightly smaller than 7 

• The decay {t + 1)"^/"^ in (7.1) can be improved easily to {t + 1)"^/^+^ for any e > 0. We 
have chosen e = 1/4 to keep the notation at a reasonable level. 

Proof. As we explained before, the proof is similar to the one in Section 3 except that now the 
function behind the front is split into a diffusive part and into an exponentially damped part 

and correspondingly there will be a few more equations. 

In Bloch space the initial conditions satisfy ||{)()IU2.« + ||^'o(' ~ ^/^' OIU^.'' < R- The 

2 , 2 ^ 

system for the variables and with initial conditions Vc\t=o = -^c^lt=0' '^slt^o — -^s^lt=0' 
and for the variable w = Wf^ ^fV with initial conditions w|i=o — ^/3,o^lt=o given in Bloch 
space by 

d,v, = + EMv^, i)J + 4K(i)„ i)J , (7.3) 

where, see (1.8) and (6.6), with v = + v^, 

M = TMT-^ , 

J^^{v^.,v^,w) = — 3C4 ■ V ®w — 3K^^ ®v ®id — V ®v ®id . 
It is useful to modify this system by introducing the coordinates (ttg, tt J by 

= ^c, «s = -^~'0U, ■v^®v^)+v^. (7.4) 

This coordinate transform takes care of the fact that asymptotically can be expressed by v^. 
Under the scaling used below the new variable Ug converges to zero, while the old variable Vg 
converges to a nontrivial expression. 

Under this transform (7.3) becomes 
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where 



We follow the lines of Section 4 and start with the renormalization process by introducing the 
scalings 

— In - 



(The 3"^ argument is the time, and the function w has here another meaning than in Section 4.) 
Note again that only the Bloch variable is rescaled, but x is left untouched. 

Under these scalings the functions ^ and still converge towards as n — > cxd. The 
variation of constant formula yields now 



+ a-'^ jW e-"'"^c,„(-r') J^^^^, v^^J) (x, X, t') , (7.6) 

— 2n 2 

+ ^-^"/^ £dr' e-""^- (X,, Jt),,,, t),, J) (x, a:, r') , (7.7) 

+ rdT'SJt,T'){^^Jv^^^,v^^^,wJ){>i,x,T') , (7.8) 
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with 

where we recall the definition 

(£/) {£, x) = f{ae, x) , 
and where S^{t, t') is now the evolution operator associated with the equation 

dj^ = a-^-{Z-M^Z-- + 7)/. . (7.9) 

Again, the exponential scaling of with respect to time does not affect the definition of 
due to the fact that only appears linearly. 

All this is quite analogous to the developments in Eqs.(4.8) and (4.9). 

7.1. The scaled linear evolution operators 

First we bound the linear evolution operators generated by M^,^ and M^^. 

Lemma 7.2. For all p e (0, 1] there exist > and 7_ > such that fori >t > r > 
and a77 (J e (0, 1) one has 

for all n e N. 

Proof. The first estimate follows directly from the fact that 

The second estimate follows from the fact that the real part of the spectrum of ^ (i) as a 
function of £ can be bounded from above by a strictly negative parabola. □ 

Next, we bound S^{t, t') as defined through (7.9) and state the analog of Lemma 4.2. 
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Lemma 7.3. Suppose that the stability condition (6.11) is satisfied. Then there is a 6 G (0, 1] 
such that for all e' e (0, 1) there exists a C^, > such that for I > r > r' > and all a e (0, 1] 
one has 

for all n e N. 

The proof of Lemma 7.3 follows closely the one of Lemma 4.2 in Section 4. 1. Therefore, 
it will be omitted here. We only remark that the estimate for the solution of (7.9) 

associated to (4.13) can be obtained exactly in the same way as (6.12). The estimates for the 
weights in £ and the derivatives with respect to x follow again as in the proof of Lemma 4.2. 



7.2. The scaled non-linear terms 

Next we estimate the scaled non-linear terms in W^^, y^^^^, and ^^.n- 

Lemma 7.4. Suppose m&x{\\v^J\^^^^,\\t\J^^^^,\\wJ\,^^J < 1. Then for all e' e (0,1) 
there exist C\ , C^, > such that for all a e (0, 1] one has 



ll>^c,nll;C^„,/, <Cl^'"'"(ll\nll;C.„ +ll^nll/C,j' 



Proof. Throughout the proof we use 

(£(/®^))(x) = a{{Zf)®{Zg))i^) . (7.11) 

i) We start with the estimates for K^^^. The most dangerous term in 

J4^{v^., v^, w) — 3U^ ■ V ®w — 3K^^ ®v®w — v®v®w 

/ 

is B^Cg^ ®v®w. From (7.11) we obtain a cr" for the scaled version ofv®w. We loose u"^ by 

taking the norm in /C^„ ^, due to the fact that is fixed and does not scale when time evolves. 

ii) We use again (7.1 1) to obtain the estimates for ^. The only difficulty stems from the 

term 

a,[M-^4(3C/, • u, ®u,)] = M-'E,{6U, ■ u, ®d^u,) 



Stability of Modulated Fronts 



32 



coming from the change of coordinates (7.4). This can be estimated in the required way by 
expressing d^u^ by the right hand side of (7.5), by using then the points ii.l)-ii.3) and the fact 
we already have a factor a"^ by ® d^u^ using again (7.1 1). 

ii. 1) The first bound for the terms on the right hand side of (7.5) is 

(with p = 1/2 for our purposes) which follows from the form of //j (£) by using the following 
lemma. 

Lemma 7.5. Let p e Cl,,i[-l/2, 1/2), C^((0, 27r), C)) with ||/i(l ■)!lc^((o,2.),c) < C\e\^^'-'^ 
forape [0, 1]. Then, there exists aC > such that for all a e (0, 1] we have 

Wi^al^Mjc^,^ < C^'^'"''^II/^IIc^,([-i/2,i/2),c2((o,2,.),c))II^IIk;, • (7-12) 



Proof. This follows since 

£2(l-p)^2(l-p) 



□ 



as 



11.2) By Lemma 7.8 below the term X^. j ^ is exponentially small in terms of a. 

11.3) From (7.1 1) we easily obtain 

iii) From [Schn96] we recall the estimates for the K^, ^ part. Note that K^.^ can be written 

^c,n = Si + S2 + ^c,n,r > 



where 

-cy^''Z-EX-''{K,u®K,u®v.,n). 
The estimate for ^NT^^ ^ follows easily by applying again (7.1 1). 
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It remains to estimate s j and §2 ■ These estimates have been obtained in [Schn96] . For com- 
pleteness we recall some of the arguments. Introducing a^(£) G Chyv^^{£,x) = a^{£)(p^n£{x) 
shows that the terms Sj and §2 are of the form 

S2ii,x) = (^a^"" Jdm Jdk K2{a''e,a''{e - m),a''{m - k),a''k) 

= (."/d™A-K.,."(.-m),."™)a„(.-™)a„(™)).„„,(.), 

with Kj : R^"*"^ — > C the kernel of an integral operator. The detailed expression for is given 
in (7.13) below. 

The case n = m — k = £ — corresponds to the spatially periodic case. In the spatially 
periodic case there exists a center manifold 

r = {u = U,^JaeR}, 

consisting of the spatially periodic fixed points related to each other by the translation invariance 
of the original Swift-Hohenberg equation. By a formal calculation it turns out that the flow of 
the one-dimensional center manifold F is determined by the ordinary differential equation 

% = • a + (0, 0, 0)a^ + ^2(0, 0, 0, 0)a^ + O(a^) . 



dt 

Since the center manifold consists of fixed points the flow a = a{t) is trivial, i.e., = 0. 
Consequently, we obtain K^(0, 0, 0) = K2{0, 0, 0, 0) = 0. Therefore, 

\K2{i,e-m,m-k,k)\ < C{\e\ + \£ - m\ + \m - k\ + \k\) , 

and so (7.11) and (7.12) imply 

Interestingly it turned out that the first derivatives of vanish as well. Since the eigenvalue 
problem M^ipg = {£)^i is self-adjoint, the projection P^{£) is orthogonal in L^(0, In) and is 
given by = {J (fg{x)u{£,x)dx)(f^{-). Thus 

K^{£,£-m,m) = 3 J dx (pg{x)(p^_^{x)(p^{x)U{x) . (7.13) 

Expanding (fii{x) = dJJ {x) -\- i£g{x) + 0{£^), with g{x) e R yields 

K^{£,£-m,m) =3 Jdx (^{d^U{x)fU{x) 



- i£g{x){dJJ{x)yU{x) + i{£ - m)g{x){dJJ{x)YU{x) 
+ {d^U{x)fimg{x)U{x) + 0{f + (£ - mf + m^) 
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Note that U{x) is an even function, so d^U is odd, which proves again (0, 0, 0) = 0. Since, 
in addition, the first order terms cancel we have 

\K^{ej-m,m)\ < C\e^ + {£ - m)^ + m^l , 

and so from (7.11) and (7.12) 

||sJ|;r < Ca^'^/^dlv llr +\\v llr )^ • 

Summing the estimates shows the assertion. □ 



7.3. Bounds on the integrals 

Here we estimate the integrals in the variation of constant formula in terms of the following 
quantities. 

Definition 7.6. For all n, we define 

K,n = sup ||i)e,n(^)llx;^„ + sup (t) ||^ , md = sup ||w^(t)||^ . 



In the following two lemmas we estimate the integrals appearing in (7.6)-(7.8). 
Lemma 7.7. Assume R^^^ + < I. Then for all\>T>a^ and all a e (0, 1] one has 

\\a-'- /Ve'^""^-(^--')(5j-,,Ji)„i)J)(-,-,r')||x:^„ < Ca^^'iRlJ , 

II /V5Jt,rO(K,,Ji)e,i)s»)(-,-,r)|L^„ < Ca-^'-^'^Rl^R:. 
Proof. We first use Lemma 7.2 and Lemma 7.4. For the second integral in (7.6) we get a bound 

sup 11.7-2- /W"^-^^-^'H3^c,n(«c,^s))(T,T) 
Te[CT2,l] Ja^ 

< Ca-'^RlJa'-/' /W'(l-r')-^/^ 



< Ca-/\RlS . 
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For the second integral in (7.7) we find similarly 



sup / V e'^""^- (3^, Ji)„i)J)(-,-,r')||;c. 



For the integral in (7.8) we find, using now Lemma 7.3 and Lemma 7.4, a bound 

< C<j<-'-^'^a^^Rl^^RZ < Ca^'-^'^Rl^^RZ . 



□ 



Lemma 7.8. Assume i?^, „ + R'Z<\. Then for alll>T>a^ and all (J e (0, 1) one has 

Proof. We restrict ourselves to the linear part Mj. A typical term of (7.6) — the first in the 
definition of Mj in (1.8) — can be rewritten as 

Note next that 

(i^ca-2"T' ®'"a-2"T')(^'^) 



= y"d£ K^^-2„^, {x-e- ip, x) w{e, X, (7-^V)e 

«N -2" ' w „^ -2-n I 



- — 2n / — 2n / 
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Using this identity, we get (because exp(Mc^^(T — r')) is bounded): 

J CT^ 

. _ — 2n / — 2n, / 

2(n 1) n 

(7.14) 

The non-Unear terms coming from can be handled in exactly the same way and yield similar 
bounds. The same is true for the terms with in (7.7). □ 

7.4. Bounds on the initial condition 

Here, we estimate the first terms on the right hand side of the variation of constant formulae 
(7.6)-(7.8). 

Lenuna 7.9. For alll >t > and alia e (0, 1] we have 
Proof. As before we have 

ll^/IL.. < ^-'/'ll/lk^„_. , (7.15) 
for < a < 1. Therefore, the first two bounds of Lemma 7.9 follow immediately from Lemma 
7.2. The third inequality is a little less obvious: First note that 

Therefore, 

||^(^„_l(Ta-^l)«;_l(■,■,l))||^^„ 

< ^-'^'\\Sn-iira-\l)w^_,{-rA)\\,c^„_, (7.16) 

< Ca-s/V-^'"e-^-""(— ')/2||^;„_i(-,-,l)||^^_^ . 

The claim is now an immediate consequence of Lemma 7.3. □ 
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7.5. A priori bounds on the non-linear problem 

This section follows closely Section 4.2. We need a priori bounds on the solution of (7.6)-(7.8). 
We (re)define now quantities analogous to those of Definition 4.3. 

Definition 7.10. For all n e N, we define 

Pcs,n = ll'^c,nlT=lllx;^„ + ll'^s,n lr=l ll/C^„ ' Pn = ll^nlT=l ll/C^n • 

Lemma 7.11. For alln there is a constant r]^> such that the following holds: If p^^ ^_ j , 
p^_i, and a > are smaller than ij^, the solutions of (7.6)-(7.8) exist for all r e [a^, 1]. 
Moreover, we have the estimates 

K,n < Ca-'pl^_, + Ce-^-""< + Ca^/\RIJ' , (7.17) 

and 

< < Ca-'/'-^'y:_, + Ca^^'-^'^Rl^K , (7-18) 
with a constant C independent of a and n. 

Remark. We remark again that there is no need for a detailed expression for rj^ since the 

existence of the solutions is guaranteed if we can show -R^^, ^ < oo and < oo. By (7. 17) and 
(7. 18) we have detailed control of these quantities in terms of the norms of the initial conditions 
and a. 

Proof. For the derivation of the estimates we assume in the sequel, without loss of generality, 
that -Rgs „ + R!^ < 1- For the first term in (7.8) we obtained in Lemma 7.9 a bound 

C(7-5/V-^'"p;;'_i . (7.19) 

For the second term in (7.8), we obtained in Lemma 7.7 a bound C(t"^^"^ ^i?^^ 

We now discuss in detail (7.7). Using Lemma 7.9 the first term is bounded by Ccr p^^ n—i- 
Lemma 7.7 and Lemma 7.8 yield for the second and third terms a bound Ca'^^'^{R^^ ^)^ + 

-n 

Ce"^'" R"^ for a C > independent of a G (0, 1] and n e N. 

Finally, we come to the bounds for (7.6). Using Lemma 7.9 the first term is bounded 
by Ca~^^^p^^ Lemma 7.7 and Lemma 7.8 yield for the second and third terms a bound 

Ca^'/^iRlJ^ + Ce-^'^'^R'^ for a C> independent of a G (0, 1] and n G N. 

The proof of Lemma 7.11 now follows by applying the contraction mapping principle to 
the system consisting of (7.6), (7.7), and (7.8). 

Then for Pes p^-i and a > sufficiently small the Lipschitz constant on the right 
hand side of (7.6) to (7.8) in C([cr^, 1], /C^.™) is smaller than 1. An application of a classical 
fixed point argument completes the proof of Lemma 7.11. □ 



Stability of Modulated Fronts 



38 



7.6. The iteration process 

As in the case of the simpUfied problem, we decompose the solution v^. ^(-j •, r) for r = 1 into 

2 

a Gaussian part and a remainder. Let V'C^^) = e"'^'*' and write 

where f^(0, a;) = 0, and the amplitude is in C. We also define 11 : /C^ — > C by 

(n/)^o = n(o)/L=o- (7-20) 

Then (7.6) can be decomposed accordingly and takes the form 

An = ^n-i+n(£dr'e-"^"^c,Ji-/)(^-2n(3^^.^ + 5^^^^ ^ (7.21) 

+ ) A„_jV(^x)^,_„^(a;) - A„V(x)^,-„^(a;) . 



If we define next pj^ = ||r„||x; + ll^s nlT=i llx; then the above construction implies p^^ ^ < 



Our main estimate is now 



Proposition 7.12. There is a constant C > such that for sufficiently small a > the solution 
(■^c,n' ■^s,n' ■"^n) (7-6)-(7.8) Satisfies for all n e N; 

\An - ^n-il < C7e-^-""< + Ca-/\RIJ' , (7.23) 



Proof. We begin by bounding the difference — A^_^ using (7.21). Since / is in H'^ as a 
function of £ we obviously have 

|n/| < C||/||^^„ . (7.26) 

Thus, it suffices to bound the norm of the integral in (7.21), but this has already been done in 
the proof of Lemma 7.7 and Lemma 7.8. 

We next bound f„ in terms of f^_i, using (7.22). The first term is the one where the 
projection is crucial: For a > sufficiently small, r^-i ^ ^a"-^ ^^^h f^_i (0) = one has 

- e-""^-(i-'^')f,_i(ax,x)||. < i\K_,\\jc , (7.27) 
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as in the proof of Proposition 4.5. This leads for the first term in (7.22) to a bound (in /C^„) 

P;-i/2 . (7.28) 
The second and third term have been bounded in the proof of Lemma 7.7 and Lemma 7.8 by 

Ce-^'^'^Rl + Ca'''^{Rlf . (7.29) 

Finally, the last term 

Xjx,a;) = e'^""^-(i-')A^_^^(^^)^^_„^(^) _ A^^(^)^^_„^(^) , 
in (7.22) leads to a bound (in /C^„): 

ll^nll < Ce-^'^""<_i + Ca^'\RlJ + Ca^Rl^^ , (7.30) 

where the last term is due to /Uj(£) = — Cj£^ + 0{l^) not being exactly a parabola. For details 
see [Schn96]. Collecting the bounds, the assertion (7.24) for follows. Finally, the bounds 
on follow the in the same way as those in Lemma 7.11. The proof of Proposition 7.12 is 
complete. □ 

Proof of Theorem 7.1. As before the proof is just an induction argument, using repeatedly the 

above estimates. Again we write C for constants which can be chosen independent of a and 
n. Assume that R = sup^gj^i?^^^ < oo exists. From Lemma 7.11 we observe for cr > 
sufficiently small, 

— 5/2— ne' w 
^ ^ Pn-l ^ ^ 5/2-ne'w 

< g^"Vs,n-i+ge-^^""< (7-31) 

< Vs,n-1 + Ce-^" Pn-1 , 

with a constant C which can be chosen independent of R. Using Proposition 7.12 we find 



cs,n < Ci\AJ+pl), 
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Therefore, we can choose u > so small that forn > 9: 

Pi <3p;_i/4 + p™_i/10 + CT--Vnl, 
Pn < Pn-l/10. 

Thus, the sequence of converges geometrically to a finite limit A^. Furthermore, we find 
that lim^^^ = 0, and lim^^^ = 0. Since the quantities \A^\, p^, increase only for 
at most 9 steps the term CR in (7.31) stays less than 1/2 if we choose \ A-^\, p[, p^ = 0{a^), 
for a sufficiently large m > 0. From (7.31) the existence of a finite constant R = sup^^j^ R^^ ^ 
follows . Finally, the scaling of w^{-,-,t) implies the exponential decay of w{t). The proof of 
Theorem 7.1 is complete. □ 



References 



[AW78] D.G. Aronson, H. Weinberger.: Multidimensional nonlinear diffusion arising in population genetics. 

Adv. Math. 30 (1978), 33-76. 
[BK92] J. Bricmont, A. Kupiainen.: RenormaUzation group and the Ginzburg-Landau equation. Comm. Math. 

Phys. 150 (1992), 193-208. 
[BK94] J. Bricmont, A. Kupiainen.: Stability of moving fronts in the Ginzburg-Landau equation. Comm. 

Math. Phys. 159 (1994), 287-318. 
[CE86] R Collet, J.-R Eckmann.: The existence of dendritic fronts. Comm. Math. Phys. 107 (1986), 39-92. 
[CE87] R Collet, J.-R Eckmann.: The stability of modulated fronts. Helv. Phys. Acta 60 (1987), 969-991. 
[CE90a] R Collet, J.-R Eckmann.: Instabilities and fronts in extended systems. 1990. Princeton, Princeton 

University Press. 

[CE90b] P. Collet, J.-R Eckmann.: The time dependent amplitude equation for the Swift-Hohenberg problem. 

Comm. Math. Phys. 132 (1990), 139-153. 
[CEE92] P. Collet, J.-R Eckmann, H. Epstein.: Diffusive repair for the Ginsburg-Landau equation. Helv. Phys. 

Acta 65 (1992), 56-92. 

[DL83] G. Dee, J. S. Langer.: Propagating pattern selection. Phys. Rev Lett. 50 (1983), 383-386. 
[Eck65] W. Eckhaus.: Studies in nonlinear stability theory. Springer tracts in Nat. Phil. Vol. 6, 1965. 
[EW91] J.-P. Eckmann, C.E. Wayne.: Propagating fronts and the center manifold theorem. Comm. Math. Phys 
136 (1991), 285-307. 

[EW94] J.-P. Eckmann, C.E. Wayne.: The non-linear stability of front solutions for parabolic partial differen- 
tial equations. Comm. Math. Phys. 161 (1994), 323-334. 

[EWW97] J.-P. Eckmann, C.E. Wayne, P. Wittwer.: Geometric stability analysis of periodic solutions of the 
Swift-Hohenberg equation. Comm. Math. Phys. 190 (1997), 173-211. 

[Ga94] T. Gallay.: Local stability of critical fronts in nonUnear parabolic partial differential equations. Non- 
linearity 1 (1994), 1A\-1(A. 

[HS99] M. Haragus, G. Schneider.: Bifurcating fronts for the Taylor-Couette problem in infinite cylinders. 
Zeitschrift jur Angewandte Mathematik und Physik (ZAMP) 50 (1999), 120-151. 

[KSM92] P. Kirrmann, G. Schneider, A. Mielke: The validity of modulation equations for extended systems 
with cubic nonlinearities. Proceedings of the Royal Society of Edinburgh 122A (1992), 85-91. 

[RS72] M. Reed, B. Simon.: Methods of Modern Mathematical Physics I-IV. New York, Academic Press, 
1972. 

[Sa77] D.H. Sattinger.: Weighted norms for the stability of travelling waves. J. Diff. Eqns. 25 (1977), 
130-144. 

[Schn94] G. Schneider.: Error estimates for the Ginzburg-Landau approximation. J. Appl. Math. Physics 45 
(1994), 433^57. 



Stability of Modulated Fronts 



41 



[Schii96] G. Schneider.: Diffusive stability of spatial periodic solutions of the Swift-Hohenberg equation. 
Comm. Math. Phys. 178 (1996), 679-702. 

[Schn98] G. Schneider.: Nonlinear stability of Taylor-vortices in infinite cylinders. Archive for Rational Me- 
chanics and Analysis 144 (1998), 121-200. 

[Ta97] M. E. Taylor.: Partial Differential Equations I: Basic Theory.Appl. Math. Sciences 115, Springer 
1997. 

[vH91] A. van Harten.: On the validity of Ginzburg-Landau's equation. /. Nonlinear Science 1 (1991), 
397^22. 

[Wa97] C.E. Wayne.: Invariant manifolds for parabolic partial differential equations on unbounded domains. 
Arch. Rat. Mech. Anal. 138 (1997), 279-306. 



